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ABSTRACT 

We present radio continuum and polarization images of the North Polar Spur (NPS) from the 
Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical 
Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency 
channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. 
Combining this RM map with a published Faraday depth map of the entire Galaxy in this direction, 
we derive the Faraday depth introduced by the NPS and the Galactic interstellar medium (ISM) in 
front of and behind the NPS. The Faraday depth contributed by the NPS is close to zero, indicating 
that the NPS is an emitting only feature. The Faraday depth caused by the ISM in front of the 
NPS is consistent with zero at 5 > 50°, implying that this part of the NPS is local at a distance of 
approximately several hundred parsecs. The Faraday depth contributed by the ISM behind the NPS 
gradually increases with Galactic latitude up to 6 = 44°, and decreases at higher Galactic latitudes. 

This implies that either the part of the NPS at 6 < 44° is distant or the NPS is local but there is a sign 
change of the large-scale magnetic field. If the NPS is local, there is then no evidence for a large-scale 
anti-symmetry pattern in the Faraday depth of the Milky Way. The Faraday depth introduced by the 
ISM behind the NPS at latitudes b > 50° can be explained by including a coherent vertical magnetic 
field. 

Subject headings: ISM: magnetic fields — ISM: structure — ISM: individual (North Polar Spur) — 
polarization — Galaxy: center — Galaxy: structure 


1. INTRODUGTION 

The North Polar Spur (NPS) is one of the largest co¬ 
herent structures in the radio sky, projecting from the 
Galactic plane at Galactic longitude / Ri 20° and extend¬ 
ing to a very high Galactic latitude b ss -|-80°. It was first 
identified in low frequency ra di o surveys in the 195 0s (e.g. 
iHanbnrv Brown et al.l llOOTil) . iLarge 01711 (jl966[) fitted 
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the NPS to part of a hypothetical circular structure with 
a di ameter of about 110° wh ich was later named Loop I 
(e.g. iBerkhuiisen et al.]ll97Hl . 

Observations and theoretical modeling of th e NPS 
up to the 1980s were thoroughly reviewed by iSalted 
(|1983f) . The NPS had by then been known to have (a) 
strong synchrotron emission whose fractional polariza¬ 
tion is very high- up t o ^70% at 1.4 GHz, at high lat¬ 
itudes (iSpoelstr al 119721) (b) to have strong X-ray emis¬ 
sion le.g. iBunner et al1ll972tl . (c) to be probably associ¬ 
ated with a vertic al Hi filament at ( ^ 40° at velocities 
aroun d 0 km s“^ (|Heiles &: Jenkins! 1 19761: iGolomb et al.l 
1980I1. and ('dl to align with starlight polarization (e.g. 


Axon &: Ellislll976ll . All these suggested that the NPS is 


an old local supernova remnant (SNR) at a distance of 
about 100 pc t hat has been reheated by the shock from 
a second SNR (jSalteilll98,‘Hl . 

There have been more observations of the NPS at vari¬ 
ous wavelengths since the 1980s. From several all-sky ra¬ 
dio continuum surveys, the brightness temperature spec¬ 
tral index {Ti, oc with Ti, being the brightness tem¬ 
perature at frequency ly) of the NPS is /3 « —2.5 between 
22 MHz and 408 MHz (iRoger et al.lll999fl an d between 
45 MHz and 408 MHz (iGuzman et al.ll2011h. and 0 w 
—3.1 between 408 MHz and 1420 MHz ( Reich fc lleichl 
[19^ at 6 > 30° where there is little contamination of 
the diffuse emission from the Galactic plane, confirming 
the NPS as a nonthermal structure. 

The NPS can be clearly seen in the soft X-ray back¬ 
ground maps from R OSAT observations, p articularly 
in the 0.75 keV band ([Snowden et al.l [199711 . Towards 
several positions, sp ectra were extracted from observa¬ 
tions with ROSAT (jEgger fc Aschenbacblll9^ . XMM- 
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Newton (|Willingale et al.ll2003f) and Suzaku ([Miller et al.l 
1200811 and fit to multiple emission components includ¬ 
ing thermal emission from the NPS. The consensus of 
these papers is that the fraction of the total Galactic 
Hi column density in front of the NPS is close to 1 
for b ^ 20° and 0.5 for b > 30°. Based on the lo¬ 
cal 3D interstellar medium (ISM) distribution from in¬ 
version o£abou^_23_000stenar light reddening measure¬ 
ments (|Lallement et al.ll2014ll and the corresponding H i 
column density distribution, iPiispitarini et all 
gued that the NPS is at a distance greater than < 

On the other hand, the N PS has also b een interpreted 
as a Galactic scale feature. iSofuel (|2000ll proposed that 
the NPS traces the shock front originating from a star- 
burst in t h e Gal actic center about 1.5 x 10^ years ago. 
ISun et ail (|2014ll showed that the lower part (6 < 4°) 
of the NPS is strongly depolarized at 2.3 GHz an d thus 
beyon d the polarization horizon of about 2-3 kpc. ISofuel 
(|2015ll found the soft X-ray emission from the lower part 
follows the extinction law caused by the Aquila Rift and 
derived a lower limit of about 1 kpc for the distance to 
the NPS, although he based this estimate on the kine¬ 
matic distance to the Aquila Rift which has very large 
uncertainties. Both these re sults suggest that the NPS 
is a G alactic scale feature. iBland-Hawthorn &: Cohenl 
(|2003f) demonstrated that the NPS can be explained by 
a bipolar wind from t he Galactic center. T here have also 
been suggestions fe.g. lKataoka et al.ll2013ll that the NPS 
is associ ated with the Fer mi Bubble (fSu et al.ll2010ll . In 
contrast, IWollebenI (1200711 modelled the NPS as two in¬ 
teracting local shells that can be connected to the nearby 
Sco-Cen association. 

A conclusive way to settle the controversy of the nature 
of the NPS is to determine its distance. In this paper we 
use radio polarization data to locate the NPS along the 
line of sight. We focus on the 1.3-1 .8 GHz polarization 
measurements fro m the Galactic Magnet o-Ionic Medium 
Survey IGMIMS: iWolleben et al.ir2010af) . By comparing 
the rotation measures (RMs) of the NPS emission with 
those of extragalactic radio sources we establish the con¬ 
tribution to Faraday depth by the ISM in front of and 
behind the NPS, and so constrain its distance. The pa¬ 
per is organized as follows. In Sect. 2, we describe the 
GMIMS data and derive the RM map, then scrutinize Hi 
and optical starlight polarization data for possible infor¬ 
mation on the distance to the NPS. In Sect. 3, we confine 
the location of the NPS and discuss the implications for 
modeling of the large scale magnetic field in the Galaxy. 
We present our conclusions in Sect. 4. 


2 . RESULTS 

2.1. The GMIMS data and the RM map 

GMIMS is a survey of the entire sky with 
spectro-polarimetry at frequencies from 300 MHz to 
1.8 GHz using several telescopes in both hemispheres 
(iWolleben et al.l[2010a[l . In this paper we use the data 
observed with the Dominion Radio Astrophysical Ob¬ 
servatory 26-m Telescope covering the frequency range 
from 1280 to 1750 MHz split into 2048 channels. Prelim- 
ina ry results from GMIMS covering the NPS were shown 
by IWolleben et al.l (|2010bll . A detailed description of 
the observations and data processing will be presented 
in a forthcoming paper (Wolleben et al. in prep.). In 
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Fig. 1.— Images of the total (upper) and polarized (lower) inten¬ 
sity of the NPS averaged between 1.44 and 1.5 GHz. The images 
are in stereographic projection centered at (Z , h) = (329°, 17?5), 
the center of Loop I IjBerkhuiisen et al.lll97lT) . a position marked 
as a blue dot in each panel. The contour marks the NPS as defined 
by its polarized intensity. The resolution is 40' and the rms noise is 
about 20 mK for total intensity and 6 mK for polarized intensity. 


summary, the observations were conducted in long scans 
along the meridian with a spacing of 12 ' to ensure full 
Nyquist sampling; a basket-weaving procedure was ap¬ 
plied to the scans to form all-sky maps at each individual 
frequency. The data have been calibrated to an absolute 
level. The final data sets are frequency cubes of Stokes 
I, Q and U with an angular resolution of 40'. 

We selected a frequency range of 1.44 - 1.5 GHz con¬ 
sisting of 253 channels where there is almost no radio 
frequency interference and derived the average total in¬ 
tensity (/ 1 . 47 ) and polarized intensity (Hi. 47 ) over this 
frequency range. The resulting images are shown in 
Fig. [T] in stereographic projection with the projection 
center at {I, b) = (329°, 17?5) which is reg arded as the 
center of Loop I (jBerkhuiisen et ahl Il971fl . We mark 
a contour denoting the outer boundary of the NPS on 
the basis of its morphology as seen in the Pi .47 image 
where Pi .47 > 0.1 K and RM errors are less than about 
5 rad m“^, as discussed below. The NPS can be clearly 
identified in both total intensity and polarization. At lat¬ 
itudes higher than about 40°, the inner edge of the NPS 
is much brighter than the outer edge, which is consistent 
with previous observations. 

For each pixel with a polarized intensity Pi .47 larger 
than 0.1 K (about her-level per frequency channel), we 
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linearly fit polarization angles (y) versus wavelength 
squared (A^) over the entire frequency range from 1280 
to 1750 MHz to obtain an RM as 

x(A')=Xo + RMA2, (1) 

where xo is a constant. The map of RMs is shown in 
Fig- m (top panel). We a lso show the Galactic Fa raday 
depth map constructed bv lOpnermann et'all (|2015h from 
RMs of extragalactic sources in Fig. [2] (lower panel). Al¬ 
though the linear fitting can also be applied to weaker 
polarized intensities, larger errors will be introduced, as 
shown below. 

It has been demonstrated that the RM from linearly 
fitting polarization angle versus A^ can be wrong unless 
the beh avior of polarization fraction again st A^ is ex¬ 
amined (jBurnlll9^ iFarnsworth et al.ll20lHl . In reality, 
the NPS is either Faraday thin with only synchrotron- 
emitting medium or Faraday thick with a mixture of 
synchrotron-emitting and Faraday-rotating medium. For 
the Faraday thin case, the linear relation between polar¬ 
ization angle and A^ always holds. For the Faraday thick 
case, the linear relation holds for certain ranges of wave- 
lengt hs and the RMs re present half of the true values 
le.g. iSokoloff et abl 1199811 . For the current data, linear 
relations between polarization angle and A^ can be seen 
for virtually all the pixels with Pi ,47 larger than 0.1 K. 
The RMs shown in Fig. [5] (top) are thus reliable. 

We also generated an RM map using all data over 
the entire frequency r ange of 1280 to 1750 MHz vi a the 
RM synthesis method (jBrentiens fc de Bruvnll200'^ . Al¬ 
though the resolution in RM is only 150 rad m“^, the 
signal-to-noise is high, allowing measurements of peak 
RM on finer scales. The resulting RM map is completely 
consistent with the RM map shown in Fig. [2] (top); the 
difference between RMs calculated in these two ways has 
a mean value of 0 ± 1 rad m“^ over an area much larger 
than the NPS. We conclude that the NPS is Faraday thin 
as the RM synthesis method ofte n fails to reprodu ce RM 
when a source is Faraday thick (jSun et al.ll2015h . The 
R M map in Fig. [2] (top) is very similar to that obtained 
by IWolleben et al.l (j2010b[l from RM synthesis based on 
the pilot GMIMS data, but has higher resolution and 
sensitivity. 

The be s t pub lished RMs for the NPS are those of 
iSnoelstr^ (jl984H . based on surveys at 408, 465, 610, 
820, and 1411 MHz. Polarization angle measurements, 
with different beams at these frequencies, were used to 
compute RM point-by-point; spatial undersampling pre¬ 
cluded smoothing to a common beamwidth and com¬ 
putation of an RM map. Because of undersampling 
in frequency, |RM| was restricted to values less than 
35 rad m“ . The resulting RM “map” is probably sensi¬ 
tive only to spatial scales >3°. No useful comparison of 
our new RMs with these older data is possible. 

We made simulations to quantify the RM errors. We 
extracted a data cube centered at an empty area with 
a size ^ 30° x 20° which contains primarily noise but 
no polarized structures in any of the frequency channels. 
For each pixel, we generated a fake source with a ran¬ 
domly selected intrinsic polarization angle, polarized in¬ 
tensity and RM, and added the corresponding U and Q 
of the source to each individual frequency channel. We 
then derived a new data cube of polarization angle and 
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Fig. 2.— RM map for the NPS region from GMIMS (up- 
per panel) and the Gal actic Faraday depth map constructed by 
lOopermann et aTI II201FI 1 from RMs of extragalactic sources (lower 
panel). The contours are the same as in Fig. [T] 
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Fig. 3. — Difference between the input and fitted RM values, 
ARM, versus the polarized intensity averaged over 1.44 - 1.5 GHz, 
Pi. 47 , from faked sources. The red shaded area shows the expected 
difference at given Pi, 47 . 

applied the same linear fitting procedure as above to cal¬ 
culate a RM map. The difference, ARM, between the 
derived RM values and the input RM values provides a 
robust estimate of the RM errors. We show ARM versus 
Pi ,47 in Fig. m We repeated the process by adding the 
same fake sources to Gaussian noise with an rms value 
of 20 mK as measured from the data, which yielded the 
expected errors (red shaded area in Fig.|3l). The real RM 
errors are much larger than the expected errors. This is 
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Fig. 4.— Hi intensity maps from the LAB survey from velocity —1.03 km to +7.21 km s“^ with a step of 1.03 km (from top 
left to bottom right). The gray scale runs from 0.1 K (white) to 20 K (dark). The NFS is outlined in red solid lines and an Hi feature 
possibly associated with the NFS is outlined in blue dashed lines. The red dots indicate the projection center as in Fig. [T] 


probably related to low-level scanning elfects in the data; 
the basket-weaving process reduces such effects, but does 
not completely remove them. We kept only those pixels 
with a Pi .47 larger than 0.1 K so that the RM errors are 
lower than about 5 rad 

Two patches with high positive values can be identified 
in both the RM and Faraday depth maps in Fig. [21 The 
one at 6 < 35° can be clearly seen in the RM map (up¬ 
per panel in Fig. [2]), but is not especially obvious in the 
Faraday depth map (lower panel in Fig. [2]). In contrast, 
the other at 36° < b < 46° is clearly seen in the Fara- 
day depth map , but ha s smaller extent in the RM map. 
IWolleben et al.l (|2010bH attribute both patches to Fara¬ 
day rotation by a local H i bubble associated with Upper 
Scorpius, one of the three subgroups of the Sco-Cen OB 
association, at a distance of about 145 pc. Towards lati¬ 
tudes above about 50°, which are not influenced by that 
Hi bubble, RMs gradually decrease to zero with large 
fluctuations. 

2.2. Hi data revisited 

iHeiles fc .lenkinj ()1976f ) claimed that the NFS is asso¬ 
ciated with Hi gas over the velocity range —20 km s“^ 


to 4-20 km s“^, using data from the Hat Creek Hi sur- 
vey. Using the Leid en/Argentine/Bonn (LAB) survey 
(iKalberla et al.ll2005[) . which has much higher sensitiv¬ 
ity, we re-investigate the associations between Hi gas 
and the NFS. 

By comparing the NFS with each individual velocity 
channel map from the LAB survey, we find that a fil¬ 
ament oriented almost parallel to I sa 40° extending 
from b ~ 30° to b 70° over the velocity range from 
— 1.03 to -1-7.21 km is possibly morphologically asso¬ 
ciated with the NFS (Fi g . |^, consistent with the find¬ 
ing by IHeiles fc Jenkinj (jl^6f ) . The vertical H i fila¬ 
ment can be best seen at velocity +2.06 km s“^, roughly 
parallel to the NFS, gradually fading away towards ve¬ 
locity +7.21 km s“^ and becoming brighter but mixed 
with large-scale background emission towards velocity 
—1.03 km s“^. A contour based on the morphology of the 
filament at velocity +2.06 km s“^ is shown in each ve¬ 
locity frame in Fig. 01 Because of the high latitude and 
the very low velocity, the distance to the Hi structure 
cannot be constrained. 

We estimate the mass of the Hi gas contained in the 
region within the dashed blue contour of Fig. [4] over the 





































The North Polar Spur (Loop I) 


5 


velocity range from —1.03 to +7.21 km to be about 
lO^H^QQ Mq, where Hioo is the distance to the Hi with 
a nominal value of 100 parse cs. Assuming t he Hi gas is 
part of a large shell structure, iWeaved (jl979ll obtained an 
expansion velocity of 2 km s“^ which corresponds to a ki¬ 
netic energy of about 4 x ergsF^ For Hioo = 1, 

the kinetic energy is 4 x 10“^^ ergs, which can be ea sily pro- 
duced by stellar winds from the Sco-Cen cluster (iWeaverl 
[T^ . and for Hioo = 100, the kinetic energ y of 4 x 10^® 
ergs is well below w hat a nuclear explosion (ISofuell2000ll 
and galactic winds l)Bland- Hawthorn fc CohenI l2003fl 
the Galactic center can provide. Thus the Hi filament 
can be either local or far away according to the energy 
budget, and the NPS, if it is associated with the Hi fil¬ 
ament, can be either local or distant. 

2.3. Optical starlight polarization revisited 

The light from stars becomes polarized when it is se¬ 
lectively absorbed during p ropagation by dust grain s 
aligned by the magnetic field (iDavis fc Greenst^ll951h . 
Starlight polarization vectors are parallel to the magnetic 
field in the dust and the polarization fraction depends 
on the depth of the sightline and on the degree of or- 
deri ng of the magnetic f ield perpendicular to the sight¬ 
line (|Fosalba et al.l[200^ . In contrast, radio polarization 
vectors, after correction for Faraday rotation, are per- 
p endicular to mag netic field vectors. 

iSooelstr^ (|1972[ 1 compared the polarization angles of 
radio emission at 1415 MHz from the NPS with those 
of optical starlight polarization and found that for stars 
with distances larger than about 100 pc the two angles 
differ by about 90° indicating that they trace the same 
magnetic field. This set the distance to the NPS at about 
100 pc. 



30 ° 0 ° 330 “ 

Galactic longitude 


Fig. 5.— Optical starlight polarizatio n vectors (red bars) overlaid 
on th e Planck dust image at 353 GHz UPlanck Collaboration et al.l 
1201411 . The lengths of the bars are proportional to the polarization 
fraction, and their orientations indicate the orientations of mag¬ 
netic fields. The green line and blue line outline the NPS and the 
Hi structure described in Sect. 12.21 The two circles are centered 
at (/, h) = (335°, 10°) with radii of 35° and 60°. 

There are now mor e optical polar i zation data in c luding 
the compilations by iHeilesI (|2000ll . ISantos et al.l (|2011il 

We cannot derive an expansion velocity from the data in Fig.[4| 
because we cannot relate this H i feature to other H i filaments to 
form a large shell structure. 
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Fig. 6. — Polarization percentage of optical starlight and the 
difference of polarization angles between the optical starlight po¬ 
larization and WMAP 23 GHz polarization, both as a function of 
distances. The filled circles (red) are for stars towards the NPS and 
the open circles (black) are for stars outside the NPS. The dashed 
lines mark a 90° polarization angle difference. 


and iBerdvugin et al.l l|2014[) . which motivate us to re¬ 
examine the correlations between starlight polarization 
and other tracers of the NPS. In Fig. [5] we show optical 
starli ght polarization data overlaid on the Planck dust 
map (IPlanck Collaboration et al.ll2014[l . 

At latitudes b > 20° the starlight polarization vectors 
have a curvature which resembles that of the dust struc¬ 
tures (Fig. [5]); the curvature of the vectors suggests a cen¬ 
ter at {I, b) Ki (335°, 10°). In Fig.[5l we show two partial 
circles with radii of 35° and 60° centered at this position. 
The starlight polarization vectors are in good alignment 
with the circles. There appears to be a dust bubble cen¬ 
tered at about the same position with a radius of about 
30°, but no prominent filamentary structure within this 
d ust bubble. _ 

iBerkhuiisen et al.l (119711) placed the center of Loop I at 
{I, b) = (329°, 17.5°), not far from the center of starlight 
polarization vectors. It is thus possible that both the 
NPS and the starlight polarization are products of the 
same field configuration. The starlight polarization vec¬ 
tors are quite well aligned with the Hi feature that we 
identify in Section 12.21 and, not surprisingly, there ap¬ 
pears to be dust associated with the H i as well. 

We conclude that the starlight polarization vectors 
cannot be firmly related to the NPS on the basis of mor¬ 
phology. We turn now to evidence provided by the per¬ 
centage polarization of the starlight and the relationship 
between starlight polarization vectors and radio polar¬ 
ization vectors (which should be orthogonal if both po¬ 
larization signals are from the same magnetic field). 

The polarization percentage of the optical starlight po¬ 
larization versus distance to the stars towards and out¬ 
side the NPS for b < 40° and b > 40° is shown in Fig. [SI 
Most of the distances are from parallax measurements 
with accuracy better than 50%. Here “towards” implies 
the area within the contour denoting the NPS in Figs. [21 
and “outside” is defined as the area 10° outside the con¬ 
tour in Fig. [21 For directions towards the NPS, we also 
show the pola rization angle differe nce from the WMAP 
23 GHz data ([Bennett et al.l 1^0131) where Faraday rota¬ 
tion is very small. 
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For b > 40°, the polarization percentages for stars to¬ 
wards and outside the NFS are very similar, and they 
both start increasing at distances above 60 pc, reach 
maximum values between 200 and 300 pc and then 
slightly decrease up to a distance of 700 pc. This can 
be interpreted by a continuous distribution of dust over 
the distance range 60 - 700 pc with the magnetic field 
inside the dust gradually changing orientation as a func¬ 
tion of distance. The angle difference is roughly centered 
at 90° for distances larger than about 60 pc, although 
the scatter is large. This indicates that the NFS traces a 
similar magnetic field to the dust, and yields a very loose 
estimate of 60 - 700 pc for the distance to the NFS. It is 
also possible that the NFS is further with its magnetic 
field extending from or coincident with the magnetic field 
in the distance of 60 - 700 pc. 

For b < 40°, there are not many polarization measure¬ 
ments for stars towards the NFS. Therefore even a rough 
estimate of the distance to the NFS is very uncertain, and 
more data are needed. 


3. DISCUSSION 
3.1. The location of the NFS 

Foll owing iBur a (fi^ and iBrentiens fc de BruvnI 
(1200^1 . we introduce the Faraday depth (FD) as a func¬ 
tion of distance along a line of sight, I, which is defined 
as ^ 

FD{1) = kJ^ ne(s)S||(s)ds, (2) 

where the integral is along the line of sight, iF is a con¬ 
stant, rie is the electron density, i?|| is the magnetic field 
projected along the line of sight, and s is the distance in¬ 
crement. For the Faraday depth of the Galaxy (FDq), I is 
the distance from the observer to the edge of the Galaxy. 
The differential Faraday depth of a source, AFD, can 
then be defined as 


AFD = FD(l2)-FD(li), 


(3) 


where h and I 2 are the distance of the near and far 
boundaries of the source, respectively. A detailed dis- 
cus sion of the distin ction between RM and FD is given 
by iSun et ahl (|2015ll . Throughout the paper, we use 
the FD map of the Galaxy which has been constructed 
bv lOppermann et al.l (1201511 pri marily based on the RM 


the Galactic Doles bv |Mao et al. 

ll20^. The extended 

critical filter (jODoermann et al.l 

2011) was used for the 


construction, which was able to simultaneously recover 
the Faraday depth, its angular power spectrum, and the 
noise co-variance. The minimum scale of the final Fara¬ 
day depth map can be as small as 0.5°. The typical 
uncertainty is about 5 rad towards latitudes greater 
than about 40°, and about 10 rad m“^ towards latitudes 
between 20° and 40°. Note that we also tried the FD 
map of the Galaxy by IXu fc Haiil (|2014ll for the anal- 
ysis below and fo u nd si milar results. The FD map by 
iQppermann et al] (j2015[l covering the NFS is presented 
in the lower panel of Fig. [2] 

Our aim is to constrain the location of the NFS by 
comparing the differential FD of the Galactic ISM in 
front of the NFS (AFDfg) with that of the Galactic 
ISM behind the NFS through to the edge of the Galaxy 


(AFDbg)- The known quantities from observations are 
the RMs of the NFS (RMnps); the FDs of the Galaxy 
through (FDg,t) and outside (FDg,o) the NFS; and 
the unknowns are AFDfg : AFDbg ; the differential FD 
of the NFS itself (AFDnps) and the differential FD of 
the Hi bubble (AFDhi). The area of the NFS has been 
outlined in Fig. [2] We restrict the area outside the NFS 
to be within 10° longitude from both sides of the NFS 
at each latitude. For Galactic latitudes between 28° and 
76° we average RMnps, FDg,t, and FDg,o in latitude 
bins of 4° and over all the corresponding longitudes and 
obtain their latitude profiles (Fig. [T]). 

The high positive RMs and FDs that are associ- 
ated with the local H i bubble in front of the NFS 
l|Wolleben et al.ir2010bl and their Figure 3) can be clearly 
seen in Fig.[2j Because of the influence of this H i bubble, 
we divided the NFS into two regions: 

• b < 50° - The differential FD of the H i bubble in 
front of the NFS has to be accounted for. We can 
represent RMnps, FDg,t and FDg, o as 

RMnps = AFDhi, t + AFDfg + ^AFDnps 
< FDg,t = AFDhi, T + AFDfg + AFDnps + AFDbg 
^ FDg, o = AFDhi, o + AFDfg + AFDbg- 

(4) 

Here AFDhi, t and AFDhi, o are the differential 
FD of the H i bubble covering the NFS and the area 
outside the NFS, respectively. The factor 4 comes 
from the assumption that the thermal gas within 
the NFS is u niformly mixed with non-thermal emit¬ 
ting gas fe.g. lSokoloff et allll998ll . The assumption 
is reasonable as good linear relations between po¬ 
larization angles and hold towards the NFS. 

• b > 50° - There is no influence by the Hi bubble 
in this region, and RMnps, FDg, t and FDg, O can 
be expressed as 

f RMnps = AFDfg + ^AFDnps 
< FDg, T = AFDfg + AFDnps + AFDbg (5) 

[ FDg,o = AFDfg + AFDbg- 

We first estimate the differential FD of the NFS. For 
the area b > 50°, it can be derived from Eq. dSD as 

AFDnps = FDg, T - FDg, o- (6) 

The results are shown in Fig. [3 (bottom panel). The 
average of AFDnps is 2±4 rad m“^, consistent with zero. 
This is likely due to the lack of thermal electrons as no 
excess Ha emission can be seen towards the NFS from 
the composite all-sky Ha map of iFinkbeineil l)2003ll . For 
the area b < 50°, the differential FD of the NFS can be 
derived from Eq. m as 

AEDnps = FDg, T — FDg,o — (AFDhi,t — AFDhi, o)- 

(7) 

The differential FDs of the Hi bubble through and out¬ 
side the NFS are unknown, it is therefore difficult to solve 
for AFDnps- For the lower part b < 36°, FDg, t— FDg, O 
is around zero (Fig. [3 bottom panel), which means 

AFDnps = AFDhi, o — AFDhi,t- 


( 8 ) 
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Fig. 7. — Galactic latitude profiles of RM, FD and differential 
FD. Top panel: the filled (open) circles represent FDs of the Galaxy 
through (outside) the NFS and the squares represent RMs of the 
NFS derived from the GMIMS data; Bottom panel: the difference 
between the FDs through and outside the NFS. 


There is no physical connection between the NPS and 
the Hi bubble and hence no relation between AFDnps 
and FDhi o ~ FDhi t- This implies that both sides of 
Eq. ® are equal to zero for the latitude range b < 36°. 
Since AFDnps is close to zero towards both larger and 
smaller Galactic latitudes, we assume it is also close to 
zero towards the middle part 36 < 6 < 50°. This can be 
corroborated by the fact that the patch with high posi¬ 
tive RM crosses the eastern edge of the NPS without any 
change (Fig. [2]). The large values of FDq, t — FDq, o for 
36 ^ ^ 50° (Fig. [7]) can be attributed to the large dif¬ 

ference between AFDhi.t and AFDhi.o which can also 
be seen from Fig. [2] For the discussions below, we assume 
AFDnps is zero for b < 50°. 

We now look at AFDbg and AFDpc. For the entire 
latitude range, we can obtain an estimate of AFDbg 
from Eqs. o and dS]) as, 

AEDbg = EDg,t — RMnps — 2 ^Fi^NPS- (9) 

From previous discussions, AFDnps is zero, and AFDbg 
can then be derived. The resulting profile of AFDbg 
is shown in Fig. [S] We can only solve AFDpG for the 
latitude range b > 50° from Eq. ([5]) as 

AFDpG = RMnps — 2^hDNPS7 (10) 

and the result is shown in Fig. [8l 

For the area b > 50°, AFDfg is consistent with 
0 rad m“^. Since AFDbg is not zero, a regular mag¬ 
netic field and constant thermal electron density must 
exist along the entire line of sight. In this case, AFDfg 
values of zero imply that the path length in the inte¬ 
gral in Eq. m is close to zero and the NPS is thus 
a local feature. From t he 3D modeling of the ISM by 
iPuspitarini et al.l (|20I4f) . the local cavity, defined as an 
irregularly shaped area of very low density gas surround¬ 
ing the Sun, extends to about 100 pc towards the NPS. 
From Fig. [51 dust starts to appear only from the dis¬ 
tance further than about 60 pc, which supports the ex¬ 
istence of the local cavity. Within the cavity, the differ- 



Fig. 8.— FDfg S'Ud FDbg versus Galactic latitude. The lowe r 
soli d and dashed lines are fro m the models by ISun Ik Reicbl 1120101 1 
and IJansson fc Farrar! I|2012|). res pective ly. The upper solid line is 
from the model bv ISun &: Reichl 11201011 but with an extra dipole 
field with a strength of 0.2 fiG at the Sun’s position. 


ential FD must be around zero, and any contributions to 
Faraday depths start beyond the cavity wall. Since the 
Galactic m agnetic field is p redominately parallel to the 
plane (e.s. ISun et al.l[20081 and references therein), the 
line of sight component of magnetic field towards lati¬ 
tudes higher than 50° is very low. The contributions to 
Faraday depths thus increase very slowly as a function 
of distance. All these considerations place the NPS at a 
distance of several hundred parsecs. 

Towards latitudes b < 50°, AFDbg increases with 
latitude from a value of —11 rad m~^ at 6 = 28° to 
17 rad m“^ at 6 = 44°, and decreases towards higher 
latitudes. There are two possible explanations for the 
behavior of AFDbg- One is that the NPS is local and 
AFDbg is from the large-scale Galactic magnetic field 
which has a change of sign at 6 = 44°. The other is 
that the low latitude part and the high latitude part of 
the NPS are separate structures. It can be seen from 
the total intensity image in Fig. [T] that the low latitude 
part is much brighter than the high latitude part and the 
transition is not smooth, which can also be seen from 
the recent all-sky radio continuum map at 1.4 GHz by 
IGalabretta et al.[(|2014f) . The comparison of polarization 
observations at 2.3 and 4.8 GHz indicates tha t the very 
low l atitude part is further than 2-3 kpc (|Sun et al.l 
l2014f) . With the low latitude part of the NPS far away, 
the path length from the NPS to the edge of the Galaxy 
is largely reduced and thus AFDbg is much less than 
that extrapolated from high latitudes. 


3.2. Modeling of the Galactic magnetic field 

Modeling of the large-scale magnetic field in the 
Galaxy is very challenging. Ideally models should be 
able to reproduce a broad range of observations such as 
the FD of the Galaxy, including the total intensity and 
polarized intensity f rom t he synchrotron em i ssion. Both 
ISun fc Reichl (1201011 and iJansson fc Farrad (|2012ll have 
built models of the Galactic magnetic field, including 
both disk and halo components. 

The differential FD of the Galactic ISM behind the 
NPS for b > 50° is almost equal to the FD of the G alaxy , 
whi ch allows us to test the rn odels of iSun fc Reichl (|2010[ 1 
and IJansson fc Farrail (|2012f) . In Fig. [51 we show the FD 
profile of the Galaxy from both these Galactic magnetic 
field models. Both models predict a FD smaller than 
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AFDbg- To increase FDs, the models need to have a 
larger magnetic field along line of sight, which can be 
achieved by increasing either magnetic field parallel to 
the Galactic plane or magnetic field perpendicular to the 
Galactic plane. From Fig. HI it can be seen that AFDbg 
tends to be constant at a value around +3 rad m~^ for 
latitudes hi gher than about 60° , consistent with the value 
obtained bv lTavlor et al.l (|2009ll from NVSS RMs for area 
76° < 6 < 90°. The magnetic field parallel to the Galac¬ 
tic plane cannot contribute FD towards the north Galac¬ 
tic pole. Therefore, a vertical component must be in¬ 
corporated to ex plain the AFDrg, w hich is not included 
in the model of ISun fc Reichl (I2010D and seems insuffi- 
cient with the X-shape magnetic field in the model of 
iJansson fc Farrad (I2012D . 

To demonstrate the necessity of a vertical magnetic 
field, we tried t o add a dipol e magn etic field component 
to the model bv ISun &: Reichl (1201011 . We find that with a 
field strength of 0.2 p,G and a direction pointing from the 
north Galactic pole towards the observer at the position 
of the Sun, the revised model can now reproduce AFDbg 
for 6 > 50° (Fig.©. 

There is uncertainty in constraining large-scale mag¬ 
netic field models with FDr g at b < 5 0°. If the 
NFS is local, the mod els by iSun fc Reichl (|2010ll and 
IJansson fc Farrad (|2012ll both fail to reproduce AFDbg- 
In this case, the differential FD of the Hi bubble domi¬ 
nates the FD of the inner Galaxy, producing the mistaken 
appearance of an anti-symmetric pattern of FDs between 
the fir st and fourth Galactic quadrants. ISun fc Reichl 
(|20I0ll incorporated this anti-symmetric pattern into 
their ov erall Galactic magneti c field model. Subse¬ 
quently, iWolleben et al.l (l20I0bll highlighted that much 
of this pattern was due to the Hi bubble, which led 
IJansson fc Farrail (j2012f ) to subtract its contribution to 
FD when modeling the Galactic magnetic field. How¬ 
ever, there still remain high FDs towards the NFS around 
b = 30° in the bottom panel of their Figure I after the 
subtraction, which is not consistent with our results in 
Fig. HI Our work demonstrates that there is no evidence 
for this anti-symmetric pattern in the large-scale FD of 
the Milky Way at least around b = 30° if the NFS is 
local. 

4. GONCLUSIONS 

The North Folar Spur (NFS), one of the largest co¬ 
herent structures in the radio sky, has been known for 
more than half a century. The nature of the NFS still 
remains controversial: is it a local supernova remnant 
or a Galactic scale feature related with a starburst or a 
bipolar wind from the Galactic center? We find that it 
can be both. 

The key to understanding the nature of the NFS is its 
location in the Galaxy, and this has been the focus of 
our paper. We employed recent Hi and starlight polar¬ 
ization data and found that neither of these datasets can 
give an exact distance to the NFS, or to the dust struc¬ 
ture within the NFS perimeter. We then turned to the 


polarization data from GMIMS for a possible constraint 
of the distance to the NFS. 

GMIMS provides an unprecedented data set with 
about 2000 frequency channels at 1.3 - 1.8 GHz. Tak¬ 
ing advantage of the multi-channel data, we were able 
to obtain an RM map of the NFS by linearly fitting the 
polarization angle versus wavelength squared. Based on 
the RM map of the NFS and the FD map of the Galaxy, 
we derived the differential Faraday depth of the NFS, 
the differential Faraday depth of the Galactic interstel¬ 
lar medium in front of the NFS and the differential Fara¬ 
day depth of the Galactic interstellar medium behind the 
NFS through to the edge of the Galaxy for the Galactic 
latitude range 28° < b <76°. 

We argue that the part of the NFS atb> 50° is local at 
a distance of about several hundred parsecs because the 
differential Faraday depth of the Galactic ISM in front of 
the NFS is around zero. This part of the NFS is proba¬ 
bly embedded in a large local magnetic field bubble that 
is traced by starlight polarization. With decreasing lati¬ 
tude, differential Faraday depth behind the NFS gradu¬ 
ally increases, reaches a maximum at 6 = 44°, and then 
slowly decreases. This implies that either the NFS at 
b < 44° is far away or the NFS is local but the large- 
scale magnetic field has a sign change. If the NFS is 
local, the large positive Faraday depths are contributed 
by an Hi bubble in front of the NFS, and the large-scale 
anti-symmetric pattern in Faraday depth is then not con¬ 
tributed by a large-scale magnetic field. 

We show t h at th e Gal actic magnetic fiel d mod els by 
ISun fc Reichl (l2010f) and IJansson fc Farrad (l2012f ) can¬ 
not reproduce the differential Faraday dep th behind the 
NFS a t b > 50°. We find that the model bv iSun fc Reichl 
(j20ir)l l plus a dipole magnetic field with a direction point¬ 
ing from the north to the south Galactic pole and a 
strength of 0.2 ^G at the Sun’s position can explain the 
differential Faraday depth behind the NFS. This demon¬ 
strates that there exists a coherent large-scale vertical 
magnetic field in the Galaxy near the Sun’s position, 
which should be taken into account in future modeling 
of Galactic magnetic fields. 

The location of the NFS is uncertain because the dif¬ 
ferential Faraday depth in front of the NFS cannot be 
solved at 6 ;< 50° due to the contamination of a local H i 
bubble in front of the NFS. Future polarimetric observa¬ 
tions at lower frequencies that provide a higher resolution 
in Faraday depth are needed to properly account for the 
Faraday depth of the H i bubble. 
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